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SUMMARY 


Data on the aerodynamic excitation of first-order vibration occur- 
ring in a representative three-blade propeller having its thrust axis 
inclined to the air stream at angles of 0°, 4.55°> and- 9 .80° are 
included in this paper. For several representative conditions the aero 
dynamic excitation has been computed and compared with the measured 
values. Blade stresses also were measured to permit the evaluation of 
the blade stress resulting from a given blade aerodynamic excitation. 

It was concluded that the section aerodynamic exciting force of a 
pitched propeller may be computed accurately at low rotational speeds. 
As section velocities approach the speed of sound, the- accuracy of com- 
putation of section aerodynamic exciting force is not always so 
satisfactory. 

First-order blade vibratory stresses were computed with satisfac- 
tory accuracy from untilted-propeller loading data. A stress predic- 
tion which assumes a linear relation between first-order vibratory 
stress and the product of pitch angle and dynamic pressure and which is 
based on stresses at low rotational speeds will be conservative when 
the outer portions of the blade are in the transonic and low supersonic 
speed range . 


INTRODUCTION 


Propeller vibration arising from inclination of the thrust axis 
relative to the air stream has recently become an important problem. 

The vibration occurs because each blade element operates at a varying 
angle of attack and Mach number as it rotates and, consequently, causes 
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fluctuating lift forces which complete an excitation cycle once each revo 
lution. Modern airplane-design trends toward higher wing loading, higher 
speed, and longer range have increased the magnitude of the change of pro 
peller inflow angle from take-off with maximum wing loading to the condi- 
tion of high speed and minimum wing loading. The total angle change for 
a given airplane may he as small as 10° or as great as 30 °, and the 
resulting stresses may he excessively high even though resonance is not 
attained. The problem has been further aggravated by the need of thin 
propeller blade sections for high-speed airplanes utilizing gas-turbine 
power plants. 

Much has been written on the subject of propeller aerodynamic 
excitation resulting from eccentric loading, and numerous methods for 
relieving blade stress caused by this type of vibration have been pro- 
posed or tried. Until recently, however, even an adequate theory for 
predicting the exciting force has been lacking, and no quantitative 
stress measurements were obtained ynder closely controlled testing con- 
ditions. In view of the seriousness of the problem, a broad research 
program has been undertaken to include the following: 

(1) Tests of a tilted propeller in a uniform flow stream at sub- 
sonic and supersonic section speeds 

(2) Full-scale tests with a wing-fuselage -propeller combination to 
determine interference effects 

( 3 ) The establishment of an adequate method of calculating the 
fluctuating aerodynamic loads (reference 1) 

(4) Investigation of the effects of blade flapping. 

This paper covers part (l) of this program. The fluctuating aero- 
dynamic loads have been measured by slipstream surveys and are compared 
with calculated values for the same conditions. In addition, stress 
measurements have been made on the blades and the results are compared 
with calculated stresses. 


SYMBOLS 


b 

Cj 

Ct* 


f dc 7 

slope of lift curve per degree 

blade width, feet 
thrust coefficient (l/pn%A) 
blade-element thrust coefficient (dCi/dx) 



NACA RM L50D13 


3 


CT'^-t instantaneous blade-element thrust coefficient 


Cfji 1 


0)t(corr ) 


instantaneous blade-element thrust coefficient corrected 
to propeller rotational speed 


c l 

c *d 

D' 

f 

f o 

h 

J 

J CDt 

H 

n 

I’cat 

q 

R 

r 

T 

t 

V 

w 

w <nt 

x 


section lift coefficient 
design section lift coefficient 
propeller diameter, feet 

natural frequency of rotating propeller blade, 
cycles per second 

static natural frequency of propeller blade, 
cycles per second 

blade-section maximum thickness, feet 

» advance ratio (v/hD) 

/ 

instantaneous advance ratio 
propeller rotational speed, rpm 
propeller rotational speed, rps 

apparent instantaneous propeller rotational speed, rps 


free-stream dynamic pressure, pounds per square foo 

6 » v2 ) 


propeller tip radius, feet 

radius to a blade element, feet 

thrust, pounds 

time, seconds 

airspeed, feet per second 

resultant velocity at a blade section, feet per second 

instantaneous resultant velocity at a blade section, feet 
per second 


fraction of propeller tip radius 
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P 

P 

lo 

0O)t 


CO 


section angle of attack, degrees 

angle of thrust axis with respect to air stream (inflow 
angle ) , degree s 

blade angle, degrees 

mass density of air, slugs per 

geometric helix angle, degrees 

instantaneous geometric helix angle, degrees 
cos a r \ 

— + sin ocrp sin cot/ 

angular velocity of propeller, degrees per second 


^tan~l 


cubic foot 
ftan' 


-1 _j\ 

JtX/ 


APPARATUS 


Propeller dynamometer .- The NACA 2000 -horsepower propeller dyna- 
mometer was used for these tests with certain minor alterations neces- 
sitated by operating the propeller with the thrust axis in a pitched 
position. Pitching of the propeller axis was accomplished by inserting 
wedges of the desired angles between the base of the dynamometer support 
strut and the tunnel floor. A nose-down position was employed and the 
propeller center was retained approximately at the center line of the 
tunnel. A sketch of the dynamometer mounted in the Langley l6-foot 
high-speed tunnel is shown in figure 1. 

Wake -survey rake .- A survey rake containing 27 to tal -pressure . tubes 
was used to obtain wake data. Figure 2 shows the survey strut installed 
in the test section. The rake was supported at one end by the dyna- 
mometer nose fairing and at the other end by the tunnel wall, with 
mounting provisions for six angular locations behind the propeller. The 
mountings were so arranged that the rake tubes were alined parallel to 
the tunnel axis. 

By use of the convention that 0° is in the direction of axis incli- 
nation, or bottom center in this case, the angular strut locations 
chosen were 75 °> 105 °, 150 °, 255 °, 285 0 , and 330 ° as measured in the 
direction of rotation. It was anticipated that the maximum blade 
loading would occur near 90 ° and the minimum blade loading near 270 °. 
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The choice of strut positions was made in order to fix as nearly as 
possible the position and magnitude of the maximum and minimum blade- 
loading points. 

Propeller . - A three -blade propeller utilizing NACA 10-(3)(08)-03 
solid aluminum-alloy blades, with NACA l6-series sections, mounted in a 
three-way adjustable -blade hub was employed for most of the tests. The 
blade-form curves, design blade-angle curve, and plan form for this" 
blade sire shown in figure 3. Some tests were also made with a two-blade 
propeller employing blades of the same design mounted in a similar four- 
way hub. In addition, a few tests were made with the same blades used 
in the latter tests mounted in a semiart iculated propeller hub. Fig- 
ure 4 is a sketch of this hub. As can be seen, the two blades are rig- 
idly connected to each other but are able to rotate as a unit about an 
axis which is perpendicular to and intersects both the propeller shaft 
axis and the blade axes. The motion is unrestrained within a limit of 
about 5° each way from the center position. 

Instrumentation . - Wire strain gages were used to measure blade 
stresses, and in all cases, were mounted on the cambered surfaces of 
the blades. Vibratory stresses were measured with Baldwin SR-4 type C 
gages and steady stress was measured with a single group of type A gages 
installed on a separate • blade . 

A schematic drawing of the wiring arrangement is shown in figure 5* 
Bridge circuits were used in connecting the gages and, in order to avoid 
temperature-compensation difficulties, each bridge was made up completely 
of strain gages mounted on the blade at the point where a stress measure- 
ment was desired. The two gages used as the active arms of each bridge 
were mounted side by side, parallel to the blade axis, and centered 
about the nominal location. The two gages used as the dummy arms of 
each bridge were mounted perpendicular to the blade axis and adjacent 
to the active arm gages. The bridge outputs were transmitted through a 
slip-ring device, or "pineapple," without amplification to the sensitive 
elements of the recording galvanometer. A potentiometer of high resist- 
ance was incorporated in each circuit for balancing the bridge. In 
order to calibrate the gage installation, access was provided into each 
circuit by means of jacks so that resistances could be paralleled with 
the two active arms of the bridge. 

In order to establish the blade angular position for peak stress 
and other trace characteristics relative to propeller angular position, 
one of the galvanometer elements was utilized as ,a revolution marker. A 
special slip ring and brush combination on the pineapple, which passed 
current for several degrees only of each complete revolution of the 
propeller, was connected to this element and a voltage applied. The 
resulting trace was a solid line except for a short break occurring once 
each revolution at the angular position of the propeller corresponding 
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to current passing through the brush and slip ring. This revolution 
marker defined the propeller angular position to an accuracy of ±2°. 


TESTS 


A blade setting of 30° at the three-quarter radius was employed 
throughout the tests and constant rotational speeds were used for all 
tests. Vibratory- stress records were made for one or more tests for 
each combination of propeller configuration, rotational speed, and pitch 
angle. Wake -survey data were recorded on all tests. 

Tests were first made with the three-blade propeller with the thrust 
axis alined with the air stream to establish a vibratory response 
arising from factors other than thrust-axis inclination and to obtain 
wake thrust data on which to base- load-vibration predictions. Subse- 
quently, the propeller was tested with its thrust axis inclined to the 
air stream at angles of 4.55° and 9 « 80 °. 

Similar tests with a two-blade propeller were made. Operation of 
this propeller, however, was restricted to 4.55° pitch angle and to 
rotational speeds not over 2000 rpm because of anticipated excessive 
dynamometer vibration. 

A limited number of tests were made with a two-blade propeller 
employing the semiarticulated hub. The hub was designed for one-blade- 
propeller research and is not properly articulated to alleviate blade 
stresses most effectively. The tests were conducted to determine whether 
the propeller disk would orient itself so as to relieve the vibratory 
stresses due to pitch. 


REDUCTION OF DATA 


Tunnel velocity .- Air-stream velocity was determined by the use of 
the tunnel calibration normally used for the dynamometer in an untilted 
condition. While the projected frontal area was increased appreciably 
by tilting the dynamometer, the cross-section-area change at any one 
position in the tunnel was so small that the effect on the tunnel cali- 
bration was negligible. The advance ratios have been corrected for 
tunnel-wall interference to free-air conditions. 

Wake -survey data .- In the reduction of wake-survey pressures to 
section thrust data, the radial and angular locations of each pressure 
tube of the survey rake were determined for each strut position and 
pitch angle by assuming that the air flowing through the propeller disk 
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continued parallel to the tunnel axis. Although this assumption is 
somewhat ''incorrect over the inboard sections, where the proximity of 
the inclined spinner undoubtedly affected the air flow, it is believed 
to represent generally a good approximation to the correct flow pattern. 
All wake-survey data have been cross-faired in order to obtain data from 
the several strut positions at identical radial locations and operating 
conditions. In the determination of incremental thrust-coefficient 
values, the data used as a zero reference were the thrust coefficients 
obtained on the untilted propeller. 

Stress data .- The strain-gage output as recorded on photographic 
paper was read directly to obtain vibratory- stress values. The stress 
was taken as proportional to the trace deflection, the constant of pro- 
portionality being determined by calibrations made on each gage instal- 
lation before and after tests. In some cases, strain-gage records were 
analyzed in detail to obtain the component^ frequencies of which they 
were composed. For such analyses four -diameter enlargement s of the 
desired parts of the records were made. The 12-point Fourier series 
method was then used to determine the components. In general, only data 
from the gage located at x = 0.45 are presented herein because the 
conclusions deduced therefrom are considered representative. 

The elements of the recording galvanometer used have a natural fre- 
quency of 100 cycles per second and are damped 0.7 of critical value. 
These characteristics fix the amount of lag of the element in record in g 
the stress variation. For the present tests, this lag varied from 18° 
to 30° for first-order vibration. In the case of the element used as a 
revolution marker, however, no lag is present at the beginning of the 
short break marking the flow of current because an element at rest begins 
to move at the same instant that current is applied. 


RESULTS AND DISCUSSION 
Aerodynamic Data 


Section-^thrust-coeff icient data . - The results of section-thrust- 
coefficient measurements for the three-blade propeller at thrust-axis 
inclination angles of 4.55° and 9.80° are presented in figures 6 and 7. 
These figures represent a summary of the measured aerodynamic data from 
these tests and are obtained by conventional fairing procedure from the 
survey-rake thrust-distribution measurements. The allowance for accu- 
racy of fairing in the value of Ct 1 should be ±0.001. For convenience 
in using the data the curves of Ct* for operation with thrust axis not 
inclined have been included in all cases. 
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A few typical examples of section-thrust-coefficient change and of 
integrated blade -thrust-coefficient change around the disk are presented, 
in figures 8 and 9 , respectively. These increments are the differences 
between values for the tilted and untilted propeller. It is estimated 
that the fairing of the six points may lead to errors as large as 5° in 
determining the values of ayt for maximum and minimum thrust 
coefficient. 

Calculated thrust-coefficient change and comparison with measured 
values . - A consideration of the basic problem of blade excitation, at 
rotational frequency, caused by inclination of the thrust axis reveals 
that the fundamental aerodynamic changes may be rather simply described. 
An analytical method has been devised by which the section aerodynamic - 
load changes may be computed from two-dimensional airfoil data (refer- 
ence 1). If elemental data for operation of the particular propeller, 
untilted, are available, obtained either from pressure-distribution 
(lift) data on rotating blades or from wake-survey (thrust) data, more 
direct methods may be used in the prediction of the load variations. 
Regardless of the method used, two fundamental considerations should be 
borne in mind. First, it must be recognized that the change in blade 
loading is influenced not only by angle -of -attack changes but also by 
Mach number changes during a revolution, and that the method used to 
predict the variations must consider both effects. These points are 
discussed more fully in the appendix. Secondly, a fundamental diffi- 
culty likely to be encountered is the lack of sufficient static airfoil 
information, so that it becomes necessary to extrapolate the data. 
Extrapolation is most troublesome at section speeds near the sonic value. 
The method used in the present report to predict load variations is 
based on elemental wake-survey thrust data obtained for this propeller 
in the unpitched condition, and is described in the appendix. 

In order to compare calculated data with measured data, section- 
thrust-coefficient change and integrated blade-thrust-coefficient change 
have been calculated for the same operating conditions as occurred for 
the measured data of figures 8 and 9 , and the results are presented in 
these same figures. Some extrapolation of the basic data used for the 
computations was necessary and has been indicated by the use of broken 
lines on portions of the calculated curves of figure 8. All of the 
extrapolation is regarded as conservative with the exception of parts 
of a few of the curves at the two most inboard stations. To simplify 
the comparison, the lag and amplitude of the thrust-coefficient change 
are considered separately. 

hag of thrust-coefficient change .- Although the aerodynamic lag is 
not expected to be of importance in the design of an actual propeller 
blade, interpretation of the results of the present investigation 
requires that a discussion be included. It is apparent that, if the 
geometry alone for the pitched propeller is considered, the maximum 
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loading should occur at the 90° position and the minimum loading at the 
270 o position with the conventions used herein. Since the method used 
to obtain the calculated values is based on geometrical considerations 
and steady-state data, these peaks occur exactly at the 90° and 270° 
positions or, in other words, show no lag. The peaks of the measured 
thrust -coefficient -change curves, however, show varying amounts of lag 
behind the 90° and 270° points. Table I includes a compilation of these 
phase differences as obtained for the complete blade from figure 9* It 
will be noted that, on the average, the minimum peak lags about 16° 
(column (8)), the maximum peak somewhat more (column (9))- Also tabu- 
lated in the table are the known causes of phase difference and their 
estimated amounts. A brief discussion follows. 

The oscillatory forces affecting the lag arise from two distinct 
causes, first, fluctuations of angle of attack in combination with a 
pulsating air flow, and second, blade flapping. It may be shown by an 
application of the equations presented in reference 2 to an airfoil 
oscillating about the 50-percent -chord position, which applies in this 
case, that the amount of lag arising from the fluctuation of angle of 
attack alone can never theoretically be greater than 5.8° and for the 
present test conditions was calculated to be 5°. Calculations have been 
made which consider the combined oscillatory effects of fluctuating 
angle of attack and flexibility or blade flapping but assume that the 
blade is a simple mass concentrated at the 0.7 station and having a 
single degree of freedom. These assumptions tremendously simplify the 
problem but undoubtedly reduce the accuracy of the results. It is 
believed that the phase lags resulting from the combined calculated 
oscillatory effects (column (4) in table I) are in the proper direction 
but that their magnitudes may be inexact. 

Because the survey rake was about 10 inches behind the propeller, 
there is some lag in the measured results due to twist in the slipstream. 
This lag was calculated at not over 0.5° at any time (column (5)). 

Another cause of lag discrepancy is based on the nonaxial flow of 
air in the tunnel. A cross flow of about 0.5° and an up flow of con- 
siderably less than 0.5° exist through the test section, which leads to 
a propeller speed excitation not in phase with that due to propeller- 
axis inclination. Calculations show that the effect will appear in the 
wake as 3° to 6° of lead rather than lag (column (6)). 

There is a difference between the phase shift for the maximum and 
minimum peak measurements (columns (8) and (9)) which may be partially 
accounted for. The difference appears to be caused, as may be noted in 
figure 8, by the many instances where the value of the section thrust 
coefficient at the 75° position was markedly low, which so affected the 
fairing as to cause a greater phase lag at the maximum peak than would 
otherwise exist. The only tangible explanation appears to be that of 
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interference by the main dynamometer support strut. In the twisting 
slipstream behind the propeller the strut becomes an airfoil at a small 
angle of attack, causing a change in the inflow angle locally at the 
propeller disk. The trend of this effect is in the correct direction 
but is insufficient to cause more than a small part of the discrepancy. 

There still exists an unexplained phase difference between calcu- 
lated and experimental angular location of peak load variation of from 
7° to 13 ° (column (12)) even with the minimum peak which is considered 
to be the more accurately located of the two. Further analysis of blade 
flapping showed that if more than a single, degree of blade freedom were 
assumed and if an allowance were made for hub flexibility and struc- 
tural damping, then the difference in the lag between estimated and 
measured aerodynamic values would be reduced by about 6°. The latter 
results indicate that the agreement should improve as the assumptions 
used in the calculations are made more realistic. 

Amplitude of thrust variations .- The agreement in magnitude of the 
computed arid measured values of applied force variation, or change in 
section thrust coefficient, appears to be good though somewhat erratic 
at the low rotational speeds (low resultant section velocity). This 
agreement may be seen in figure 8 in the data for rotational speeds of 
1350 and 1600 rpm. As the computed values are based on simple geometric 
considerations and steady-state data, it may be concluded that consid- 
eration of oscillatory lift forces is unnecessary. At the higher rota- 
tional speeds (2000 and 2160 rpm) the agreement is not so good. It has 
proven difficult to ascribe the discrepancy to any particular phenomenon. 
A possible explanation for the discrepancy at the higher speeds is that 
the dynamic lift-curve slope may be appreciably different from the static 
lift -curve slope for the same average Mach number. In this event the 
computed values would be in error since they are based on steady-state 
data . 


,It seems doubtful, however, that such a difference between the 
static and dynamic lift-curve slopes would be a predominant factor in 
the discrepancy. These lift-curve-slope differences would be expected 
to be dependent on the reduced frequency, which is the flutter parameter 
commonly used to relate the frequency of oscillation to the blade chord. 
Throughout these tests the reduced frequency was essentially a constant 
and, thus, differences between static and dynamic lift-curve slopes, if 
present, should have been apparent at all rotational speeds. 

The calculated values of figure 9 for blade-thrust-coefficient 
variation show much better agreement with the measured variation than 
do the section data, especially with regard to the over-all or peak-to- 
peak values. The discrepancies of the section data tend to average out 
when integrated to obtain blade thrust. 
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To show further the agreement in magnitude between computed and 
measured values of blade thrust variation, one-half of the peak-to-peak 
displacements of the thrust-coefficient-variation curves shown in fig- 
ure 9 have been converted to pounds of thrust and plotted in figure' 10. 
The close agreement is apparent. For the abscissa in this figure, the 
product of the pitch angle and the free-stream dynamic pressure a/pq 
has been used, since it may be shown that the exciting force theoreti- 
cally varies directly with this product. The linearity is illustrated 
by the straight line drawn from the origin. A derivation of the parame- 
ter arjiq may be found in reference 3» 

It might be thought that the severe effects of compressibility 
encountered at the higher rotational speeds would destroy the linearity 
in figure 10. However, the compressibility effects, while severe, 
result in a shift of the loading toward the inboard sections and do not 
seriously affect total blade thrust variation as shown in figure 10. 

The reapportionment of the load between the inboard and outboard sec- 
tions is illustrated in figure 11 where one, -half of the measured peak- 
to-peak section-thrust -coefficient change for each blade station is 
plotted against blade radius for four rotational speeds at the same 
inclination angle. 


Stress Data 

Data from the strain gage located at the x = 0.45 station were 
selected for detailed analysis because it was the approximate location 
of maximum blade stress and because data obtained therefrom were con- 
sidered to be representative. It should be recognized that the stresses 
measured by the gage located at x = 0.4-5 represent the effect of 
vibratory forces on all of the blade outboard of x = 0.45 and. not for 
that station alone. 

An inspection of the vibratory stress records obtained in the cur- 
rent tests indicated the presence of excitations at frequencies other 
than the principal propeller-speed excitation frequency. It was neces- 
sary, therefore, to perform a harmonic analysis of the vibration records 
to isolate the once-per-revolution, or 1-P, vibration component. The 
assumption that the resulting first-order stress curve is a true repre- 
sentation of the stress variation due to inclination of the thrust axis 
cannot be fully justified because the harmonic analysis performed derives 
a first-order stress variation which is a true sinusoid, whereas the 
expressions derived for the thrust-coefficient variation (see appendix) 
indicate that this excitation is not exactly sinusoidal. 

The extraneous orders of stress excitation were largely twice -per- 
revolution, or 2-P, which was believed to be an excitation arising prin- 
cipally from the dynamometer. At the lowest rotational speeds this 
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2-P excitation closely approached resonance. Higher-order frequencies 
were noted but were generally small in magnitude. 

Lag values .- In order to obtain stress lag data which could be com- 
pared with the lag values from figure 9 as tabulated in table I , (col- 
umn (10)), the traces analyzed were from test points for which the oper- 
ating conditions duplicated as closely as possible those of figure 9- 
The lag values thus obtained from the 1-P component of the trace, cor- 
rected for galvanometer element lag, are included in table I (column (ll)). 

The phase shift between impressed force (measured aerodynamic 
average) and the analyzed 1-P trace maxima is from -1° to 12° (col- 
umn (13)) and averages about 6°. The agreement is considered reasonable 
because of the accuracy of interpretation of timing (±2°) and of aero- 
dynamic phase shift (±5°). 

Magnitude of stresses .- It was desired to investigate the trend in 
measured stresses by using the same parameter a/pq as was used for the 
exciting force in figure 10. The test conditions chosen to obtain the 
points for this plot (fig. 12(a)) approximate peak propeller efficiency 
conditions but do not duplicate those for the thrust comparison plot of 
figure 10. For some combinations of rotational speed and thrust-axis 
inclination, data from two points near peak efficiency were analyzed so 
as to indicate the shift in stress with a slight change in advance ratio. 
The stress points at the low rotational speeds conform to a linear rela- 
tionship} however, the stresses obtained for the higher rotational 
speeds (2000 and 2l6o rpm) are lower than would be indicated by a linear 
extension of the low-speed values. The lower stress values at the 
higher rotational speeds were to be expected because the trend already 
shown in figure 11 indicated an inboard movement of the change in load 
with increased rotational speed, and the stresses at the blade station 
considered would thus be reduced below that anticipated by a linear 
extrapolation. This fortunate attenuation of the stresses when t he 
outer portions of the blade are in the transonic and low- supersonic 
speed range will lead, therefore, to a conservative design if the design 
is based on the conditions existing at the lower rotational speeds. 

When all sections are operating supersonically, this stress attenuation 
will become less because the blade loading will revert to a loading 
similar to that found in the subsonic case. 

In order to establish that this stress reduction is associated 
with the inboard movement of load with increased rotational speed, a 
blade vibratory stress calculation based on the measured aerodynamic- 
load variations was made for each of the experimental values of fig- 
ure 12(a). The method of calculation involved finding by an iterative 
process the correct moment distribution along the blade. At the two 
lower rotational speeds, both measured and calculated thrust-change data 
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have been used to compute blade stresses) however, at the two higher 
rotational speeds the stress values have been determined only from 
measured values of section-thrust -loading change because of the exces- 
sive extrapolation required to -obtain calculated thrust -loading change. 


The stress values determined by- the calculations are plotted in 
figure 12(a) and indicate considerably lower stresses than were measured 
because the magnification of the stresses associated with an approach 
to resonance has not been accounted for. The natural frequency at each ' 
rotational speed was then calculated by the usual equation 
(f^ = fg + cn^) by using an assumed centrifugal-correction factor c 
of 1.7 for the first mode and a static natural frequency determined by 
bench tests to be 28.8 cycles per second. The magnification factor at 
each rotational speed was determined, at the proper ratio of propeller 
speed (l-P) frequency to natural frequency, from the conventional 
response curve for forced vibration of a system with a single degree of 
freedom and approximately zero damping. The computed stresses when cor- 
rected by this amplification factor show favorable agreement with the 
measured stresses (fig. 12(b)). Relief in the magnitude of blade thrust 
variation caused by the inboard movement of the load is verified by the 
reduction in stress at a rotational speed of 2160 rpm. 


So far only peak propeller efficiency conditions have been con- 
sidered. The parameter aq>q must be modified in any comparison of 
stresses for widely divergent operating conditions. The parameter 
should take into account the actual lift coefficient, geometric helix 
angle, and lift-curve slope. The full expression relating the lift 
variation (exciting force) to these variables has been derived 
(reference 3) I 

^inax = a T c l( a + 2c l cot 0 o ) b ^ 

The factor cu^q was therefore multiplied by (a + 2cj cot 0 o )b to 
obtain the modified form, the theoretical lift-curve slope corrected for 
Mach number being used rather than the undetermined actual values. The 
stresses plotted in figure 13 against the resulting modified a^q are 
obtained from unanalyzed data, but the error introduced by the omission 
of a harmonic analysis was investigated and found not to affect the com- 
parison of the resulting curves. The net effect of the modification 
.of ctq.q is to give a linearity to all of the stress data at the low 
rotational speeds. The comparison indicates the same trend, however, 
as for the peak efficiency points plotted in figure 12. It should be 
noted that this comparison is largely dependent on the choice of a rep- 
resentative blade station, x = 0.7 in the present case) however, no 
station appeared to give entirely satisfactory results as far as 
achieving the expected linearity of the stress - cnjq relation for all 
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speeds. The conclusion may therefore he drawn that lift coefficient 
and lift-curve slope at high section Mach numbers are not given full 
consideration in the method of application. 

Comparison of stresses for two- and three-blade propeller .- A com- 
parison is made in the following table of blade stress at x = 0.5 for 
the three configurations investigated at several rotational speeds; 
three-blade, two blades in rigid hub, and two blades in semiarticulated 
hub : 


Propeller 

rotational 

speed 

(rpm) 

Stress at x = 0.5 
(lb/sq in. ) 

3 -blade 

2 -blade 

2 -blade 

semiarticulated 

1350 

794 

840 

528 

1600 

1008 

1100 

460 

2000 

1203 

i 

1256 

— 


The stress values were all obtained at propeller -thrust -axis inclina- 
tions- of 4.55° and are for value s of J which represent the advance 
ratio for peak efficiency of the three-blade propeller. 

No harmonic analysis has been performed on the traces for which 
the stress values are presented because of the varying percentages of 
1-P stress indicated by the original records. The table compares, 
therefore, vibratory stress without regard to the order. The stresses 
for the three -blade propeller were predominantly 1-P and were of the 
same order of magnitude as those for the two -blade rigid propeller, 
which, however, contained somewhat more of the 2-P component. The semi- 
articulated two-blade propeller indicated considerably reduced stresses 
from those for a rigid propeller. The reduction indicates a reorienta- 
tion of the propeller plane in such a manner as to reduce greatly the 
first-order stresses but is accompanied by an increase in higher-order 
stresses, principally 2-P. At least a part of this 2-P vibration was 
probably caused by the type of articulation incorporated in the hub. 

With this type of articulation, an oscillatory component is superimposed 
on the steady component of rotational speed which is transmitted to the. 
propeller. This oscillation is of twice -per-revolution frequency and 
will be present whenever the propeller plane is not normal to the shaft. 

The example as cited represents an approach to stress reduction 
through thrust axis tilt, as with flexible engine mounts, rather than a 
solution by articulation. The propeller plane in this case was free to 
orient itself as much as 5° out of a plane perpendicular to the shaft. 
The propeller appeared somewhat unstable as evidenced by considerable 
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.variation in the stresses. The instability was probably caused by 
slight changes in tunnel conditions and in magnitude and point of appli- 
cation of the forces acting on the blade. The articulated blade used • 
represents an improper approach to articulation as a solution to the 
1-P stress problem. An articulation in the form of an angle -of -at tack 
change, such as is incorporated in helicopter rotors, should be of 
greater benefit. • 


CONCLUSIONS 


The following conclusions may be drawn from the present tests: 

1. The aerodynamic exciting force that results from a pitched 
attitude of the propeller may be calculated with adequate accuracy at 
low rotational speeds. 

2; The accuracy of the calculated values of section aerodynamic 
exciting force when section speeds approach sonic is not always as sat- 
isfactory as at the low rotational speeds. 

3. In the prediction of aerodynamic -load changes of a propeller 
inclined to the air stream, it is desirable to use section data in a 
manner to allow for the change of characteristics with change of Mach 
number around the disk. Such a procedure becomes necessary when sections 
operate at transonic speeds. 

4. First-order blade vibratory stresses may be computed with satis- 
factory accuracy from untilted-propeller loading data. 

5. A stress prediction which assumes a linear relation between 
first-order vibratory stress and the product of thrust-axis inclination 
and dynamic pressure, and which is based on stresses at low rotational 
speeds, will be conservative when the outer portions of the blade are 
in the transonic and low supersonic speed range. 


Langley Aeronautical Laboratory 

National Advisory Committee for Aeronautics 
Langley Air Force Base, Va. 
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APPENDIX 

A METHOD OF CALCULATING AIR LOAD VARIATIONS 
ON A TILTED PROPELLER 


General description of method . - A consideration of the geometry 
of a propeller operating with its thrust axis inclined to the air stream 
permits determination of the instantaneous values of J and other vari- 
ables fixing conditions on the blade. If the values thus obtained are 
applied to the known loading characteristic curves of the same propeller 
for the untilted condition, it is possible to predict the instantaneous - 
loading on the blade at any desired point in the cycle of load variation 
which takes place in each complete revolution. This procedure neglects 
the dynamic effects of a changing angle of attack and a pulsating air 
velocity on airfoil characteristics but, since the wave length of these 
variations will be a number of chord lengths, the effects should be 
small. It also neglects the inertia effects of the changing physical 
distortion of the blade which will be associated with forced vibration 
of the propeller. 

Velocity diagrams .- A velocity diagram showing the three-dimensional 
character of the velocity variation affecting the blade section is shown 
in figure 14(a). Induced velocities have been neglected. Relative to 
the blade section, the air-stream velocity vector V moves so as to 
cause the resultant velocity vector W to vary not only in magnitude 
but in direction. The directional variation effects a cyclic change in 
both the angle of attack and the sweep angle. 

If all velocity components parallel to the blade axis are neg- 
lected, that is, those components causing the sweep-angle variation, a 
simplified diagram may be drawn, as in figure ll+(b). In order to retain 
the customary convention which requires an air-stream velocity vector 
at right angles to the rotational velocity vector, the blade element is 
considered to be operating .with an air- stream component of V cos clt 
constant in direction and magnitude, and with a rotational compo- 
nent jmDx + V sin cup sin out constant in direction but varying in 
magnitude as the propeller rotates. 

The apparent instantaneous value of the rotational speed n^t 
may then be determined: 

jtn^Dx = iraDx + V sin a<p sin cut 
V sin oj. sin cot 


^t = n + 


nDx 
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The instantaneous value of the advance ratio J^t will be 


ItX COS Orji 


00 t rtx 

-j + sin Oiji sin o>t 


The diagram of figure 14(b) and the value of derived there- 

from are identical with those of reference 1. 


Construction of plots for the determination of instantaneous 
values of Cp'.- It is now necessary to convert the section thrust data 
for the untilted propeller into a form suitable for the determination 
of the instantaneous values. Such a form is represented by the plots of 
figure 15 , in which the section thrust-coefficient values for several 
values of rotational speed are plotted for the untilted propeller and 
the curves of constant J are cross-faired on the plots. In this form, 
the instantaneous value of section thrust coefficient Ct' at any 

desired condition within the scope of the plot and for any angular posi- 
tion of the blade may be read directly when the forward-velocity compo- 
nent V cos cup and the instantaneous advance ratio J^t are known. 


The value of C^' thus read, however, is based on an apparently 
cot 

varying rotational speed and must be corrected to the true rota- 

tional speed n. The conversion is readily made (reference 1): 


^ cot(corr. >) 



= Op*^^! + — sin ct/p sin cut^ 


Attention should be called to the effects of supercritical section 
velocities on the air load variations obtained by tilting the propeller. 
As shown on the plots of figure 15, the thrust coefficient for constant 
advance ratio drops rather sharply as the rotational speed is increased 
beyond 2000 rpm for sections outboard of the x = 0.7 station. This drop 
corresponds to the reduction of lift coefficient for constant angle of 
attack experienced with most airfoils as the speed is increased beyond 
the force -break Mach number. The effect on the section thrust- 
coefficient variation and thus on the resulting vibration is that of 
relief, in that the variation experienced is greatly lessened over that 
which would have been expected had the drop not occurred. It is there- 
fore important that the variation of lift coefficient with Mach number 
be considered as well as the lift-curve slope. There are probably 
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several methods other than the one presented herein by which these 
effects can be included. It is important, however, that both variables 
be accounted for. • The results obtained by a consideration of the lift- 
curve slope alone will be satisfactory for the lower section speeds 
where the change of lift coefficient with Mach number is small. In the 
region where the value of the lift coefficient is changing rapidly with 
Mach number, the predictions will be considerably in error. 
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Figure 1.- Eynamometer and survey strut. 
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NACA 


Figure 2.- Close-up of survey rake. 
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Figure 3.- Blade-form curves for NACA 10-(3) (08) -03 propeller. 
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Figure k.- Schematic diagram of semiarticulated propeller hub 
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(b) x = 0.^5. 

Figure 6.- Section thrust coefficient for three-blade propeller, aip = 4.55°. 
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Figure 6.- Continued 
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Figure 6.- Continued 
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(f) x = 0 . 80 . 
Figure 6.- Continued. 
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Figure 6 


Concluded. 



Section thrust coefficient , 


38 


NACA RM L50D13 




10 12 14 



(a) x = 0.30. 

Figure 7*- Section thrust coefficient for three-blade propeller, cup = 9*80° 
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Figure 7 .- Continued 
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Figure 'J.- Continued 
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Figure 7 -- Continued 
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Figure 7»- Concluded 
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Change in section thrust coefficient, AC 
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Figure 8.- Continued. . 
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Figure 8.- Continued. 
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(f) a T = 9.80°, 2000 rpm, J = 1.25. 
Figure 8.- Continued. 
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Figure 8.- Concluded. 
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Figure 10.- Thrust variation as a function of the product a^q 
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(a) Calculated values not corrected for resonance amplification. 
Figure 12.- Experimental and calculated stress comparison. 
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(i) x = 0.95* 
Figure 15 .- Concluded. 
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